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ABSTRACT 

One way of recovering information about the initial conditions of the Universe is 
by measuring features of the cosmological density field which are preserved during 
gravitational evolution and galaxy formation. In this paper we study the total number 
density of peaks in a (galaxy) point distribution smoothed with a filter, evaluating its 
usefulness as a means of inferring the shape of the initial (matter) power spectrum. 

We find that in numerical simulations which start from Gaussian initial conditions, 
the peak density follows well that predicted by the theory of Gaussian density fields, 
even on scales where the clustering is mildly non-linear. Eor smaller filter scales, 
r < 4 — 6 /i“^Mpc, we see evidence of merging as the peak density decreases with time. 

On larger scales, the peak density is independent of time. One might also expect it to 
be fairly robust with respect to variations in biasing, i.e. the way galaxies trace mass 
fluctuations. We find that this is the case when we apply various biasing prescriptions 
to the matter distribution in simulations. If the initial conditions are Gaussian, it 
is possible to use the peak density measured from the evolved field to reconstruct 
the shape of the initial power spectrum. We describe a stable method for doing this 
and apply it to several biased and unbiased non-linear simulations. We are able to 
recover the slope of the linear matter power spectrum on scales k < 0.4 /i“^Mpc“^. 

The reconstruction has the advantage of being independent of the cosmological 
parameters (D, A, Hq) and of the clustering normalisation (erg). The peak density and 
reconstructed power spectrum slope therefore promise to be powerful discriminators 
between popular cosmological scenarios. 


Subject headings: Large-scale structure of the universe — galaxies: clustering, methods: 
numerical, statistical 
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1. Introduction 

Attempts to reconstruct the initial conditions from which large-scale structure in the Universe 
grew have been carried out using several different methods. One avenue of research involves 
attempting to “run gravity backwards” from the present day galaxy distribution. This has been 
done with dynamical schemes (eg. Peebles 1989, Nusser &: Dekel 1992, Croft & Gaztahaga 1997) 
and also by applying a one to one mapping between the final smoothed density and the initial 
density (assumed to have a Gaussian p.d.f.) (Weinberg 1992). A different approach attempts to 
recover statistical measures of the initial conditions (eg. the two-point correlation function or the 
power spectrum) from knowledge of how gravitational evolution has affected these statistics. For 
example, formulae have been proposed which map the correlation function measured at the present 
day onto the initial correlation function (Hamilton et al. 1991). The same is true for the power 
spectrum (Peacock and Dodds 1994, Jain, Mo & White 1995 , Baugh & Gaztahaga 1996). In this 
paper we will show how the initial power spectrum shape can be inferred from measurements of 
the space density of peaks in the galaxy distribution, even on scales where significant non-linear 
evolution has taken place in the underlying mass. 

Peaks in the initial density distribution have been thought to be potential sites for the 
formation of galaxies and clusters of galaxies. Two seminal papers by Kaiser (1984) and Bardeen 
et al. (1986) in which the details of the theory of Gaussian random fields relevant to these problems 
were studied have been very influential in the study of galaxy and large-scale structure formation. 
In the present paper we apply these results to the density field smoothed on larger scales, in order 
to see how peak theory describes the properties of an evolved density field. We are particularly 
interested in seeing whether the distribution of peaks enables us to recover some information about 
the initial conditions, and also in its potential use as a discriminator between cosmological models. 

When applying peak theory to the problem of galaxy or cluster formation, the assumption 
is usually made that these objects have formed from the gravitational collapse of matter around 
a high peak. Some smoothing scale, somewhat arbitrarily chosen, is associated with the object 
to be formed. The validity of these assumptions has been tested on the small scales relevant to 
the formation of individual objects by Park (1991), and Katz et al. (1993) amongst others. The 
conclusions seem to be that at least for the formation of matter haloes, there is only qualitative 
agreement between the predictions and the results of numerical simulations. Indeed in view of 
the extensive merging of matter clumps predicted by presently favoured hierarchical models of 
structure formation (see e.g. Press &: Schecter 1974, Lacey &: Gole 1993) it is reasonable to expect 
that evolution of the density field on small scales will tend to disrupt the predictions of peak 
theory. 

When the matter distribution is smoothed on somewhat larger scales however, for example in 
the estimation of the genus (e.g. Melott, Weinberg & Gott 1988), the topology and some aspects 
of large scale features of the density field seem to be little affected by gravitational evolution. In 
particular, the rank order of density contrasts measured at points in the field is approximately 
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conserved. This fact has been used by Weinberg (1992) in the reconstruction method mentioned 
above. This reconstruction method also appears to work well when tested on simulations where 
the matter distribution is biased using some ad hoc galaxy formation model. 

In this paper we will use numerical simulations to test whether this insensitivity to 
gravitational evolution affects the mean number of peaks and if so on what scales it begins to 
break down. We are also interested in reconstructing the power spectrum shape from the peak 
density as one specifies the other entirely in Gaussian models. The outline of the paper is as 
follows; In Section 2 we will briefly introduce the concept of the number density of peaks in 
Gaussian random fields and in Wbody simulations. We will describe the simulations and our 
method of peak selection before testing peak theory against the evolved density both of mass and 
of ‘galaxies’ identified in the models. In Section 3 (and Appendix A) we will describe a stable 
method of recovering the power spectrum slope from the peak density. We will test it first on 
analytic power spectra and then on the simulation results. In Section 4 we discuss our results and 
present some conclusions. 


2. The number density of peaks 

We consider local maxima of all heights in a density field to be peaks. In this paper we will 
only consider cosmological models which start from Gaussian initial conditions. The derivation 
of quantities related to the power spectrum of fluctuations will rely on this assumption, but of 
course the raw peak density itself when measured in observations can be used to apply constraints 
to other sorts of models. 


2.1. Peaks in Gaussian random fields 

For a Gaussian random field the number density of peaks per unit volume npk is given by 
(Bardeen et al. 1986): 
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In a practical situation the field has to be smoothed on the scale, r/, of interest. We will use a 
Gaussian window filter, so that P{k) in the above equations should be multiplied by exp [— 


2.2. Peaks in numerical simulations 

To find maxima we need the density field to be differentiable, which means that a filter 
must be applied to the particle based realisation of density derived from Wbody simulations. An 
additional limitation imposed on us by simulations is their limited ability to resolve underdense 
regions, which will largely affect minima. These are present in identical numbers to maxima in 
Gaussian fields. The sampling of maxima being in general better than that of minima (most local 
maxima are ‘higher’ than most local minima) we choose to study only maxima here. 

We use simulations of three different spatially flat cold dark matter dominated Universes. One 
set of simulations is of “standard” CDM, with = l,h = 0.5, and another is of low density CDM 
with flo = 0.2, h = 1 and a cosmological constant A = 0.8 x 3i?o- The power spectra, P(k) for 
these models are taken from Bond & Efstathiou (1984) and Efstathiou, Bond and White (1992). 
The shape of P(k) is parametrised by the parameter P = fl/i, so that we have P = 0.5 CDM and 
P = 0.2 CDM. We also use simulations of a Mixed Dark Matter universe where a massive neutrino 
component contributes = 0.3, ^cdm = 0.6 and ^haryon = 0-1. The power spectrum was taken 
from Klypin et al. (1993). Each simulation contains 10® particles in a box of comoving side-length 
300 /i.“^Mpc and was run using a P^M Wbody code (Hockney and Eastwood 1981, Efstathiou 
et al. 1985). The mean comoving interparticle separation is therefore 3 /i.“^Mpc, of the same order 
as that of normal galaxies. The simulations are descibed in more detail in Dalton et al. (1994). 
In this paper we use 5 realisations of each model, error bars being estimated from the standard 
deviation of their results. 

To find peaks, we first assign the particle density to a 256^ grid using a cloud in cell scheme ( 
Hockney & Eastwood 1981). We then smooth the density field in Fourier space using a Gaussian 
filter of the appropriate radius. We then identify local maxima on the grid in real space. We carry 
out this procedure for several logarithmically spaced values of the filter radius. We find that if 
the filter radius is smaller than the mean interparticle separation the number of peaks increases, 
due to discreteness effects, as isolated particles become erroneously identified as peaks. The finite 
grid may also have a small effect in the opposite direction as we cannot identify peaks which are 
less than two grid cells apart. In this paper we choose to analyse fully sampled simulations, using 
filter radii from 3 /;.“^Mpc to 50 /i“^Mpc. We defer the study of discreteness effects, which will be 
important in the study of real galaxy surveys, to future work. In finding peaks, we also make use 
of the periodic boundary conditions in the simulations. The boundaries of real surveys will also 
have an effect which will need to be quantified. 

Our results, the number density of peaks as a function of filter radius are shown in Figures 
1 and 2. When plotting our results, in order to make the figures clearer, we have chosen to plot 
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Fig. 1.— Peak space density (actually Upkr^, proportional to the mean number of peaks per 
smoothing volume) for three different output times in the Standard (F = 0.5) CDM scenario. The 
errors on the points are calculated from the standard deviation about the mean of results from 5 
simulations. The linear theory prediction is shown by a curve. 



Fig. 2.— Upkrjr for three different cosmological models, F = 0.5 CDM, F = 0.2 CDM, and Mixed 
Dark Matter The linear theory predictions are shown by curves. 






- 6 - 


not the peak space density but Upk multiplied by This is proportional to the number of peaks 
enclosed on average by a smoothing volume. We show results for the T = 0.5 CDM simulation at 
different stages in its evolution (Figure 1), as well as results for the other models (Figure 2). The 
simulation points are labelled according to the predicted linear theory amplitude of fluctuations in 
8 /i“^Mpc radius spheres, cjs. The models shown in Figure 2 have a fluctuation amplitude in the 
range allowed by the 4 year COBE results (Hinshaw et al. 1996). In the figures, we also show as 
curves the predictions of peak theory of Gaussian random fields calculated using the formulae of 
Section 2.1. 

During the course of gravitational evolution, we do expect the peaks to move and to change 
in height through accretion of matter. Indeed, the number density of objects above a certain mass 
threshold changes rapidly, and can be used as a sensitive probe of the value of ug (see eg. White, 
Efstathiou and Erenk 1993). However, in this paper we are interested in the total space density of 
peaks of all heights, which is predicted to be constant in time in linear theory. 

We can see from the first two figures that the measured number of peaks follows the theoretical 
prediction remarkably well, particularily on filter scales rj > 4 — 6 /i“^Mpc. We attribute to 
merging the slight decrease in the number of peaks over time which we see in Figure 1. Of the 
two CDM models, the F = 0.5 model, which has relatively more power on small scales departs the 
most from the prediction. Merging is not a very important process, though, as even on scales of 
3 /i“^Mpc, the space density of peaks only falls ~ 10 — 15% below the peak theory prediction by 
the time the simulation has reached cjg = 1. There appears to be even less merging in the MDM 
model. 

The amplitude of fluctuations in the density field smoothed on the scales where the peak 
theory prediction begins to break down is < <5 >~ 1 — 2, which is in what is conventionally referred 
to as the quasilinear regime. We are therefore probing roughly the same scales where higher order 
perturbation theory (eg. Baugh, and Efstathiou 1994) can be used to estimate departures from 
linear theory. 


2.3. Biased simulations 

The relationship between galaxies and mass is one of the great uncertainties which affects 
the study of large-scale structure. It would be very useful to be able to measure properties of 
the galaxy density field and know exactly how they are related to those of the underlying mass. 
In the calculation of statistical clustering measures such as the two-point correlation function the 
uncertainty is parametrised with a linear biasing parameter. The value of this parameter in the 
real Universe, is still an unknown, if indeed it can be measured at all. In the case of the peak 
density, we can be hopeful that differing relationships between galaxies and the mass will have a 
smaller effect than on the two-point correlation function. It seems reasonable to expect that where 
a concentration of matter forms a peak in the density field there will also be a peak of some sort 
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Fig. 3.— The effective bias b = a galaxies (r) / amass (r) as a function of scale for 4 different biasing 
prescriptions applied to the F = 0.5 simulations (see Section 2.3). 

in the galaxy distribution. In this case, when we measure the total space density of peaks in the 
galaxy distribution we are measuring the equivalent quantity for the mass. Of course it is possible 
to imagine that no correlation exists between the galaxies and mass, or that the relationship is 
extremely complex and very non-local, so that our assumptions do not hold. Fortunately, there is 
some evidence that this is not the case, at least on large scales. For example, there appears to 
be a good correlation between galaxy velocities predicted from the density field and the actually 
measured velocities, which respond to the underlying mass (see eg. da Costa et al. 1997). 

In order to partially test these ideas, we have applied different ad hoc biasing prescriptions 
to the F = 0.5 simulation in order to produce a set of galaxy particles. One prescription involves 
carrying out the following procedure. For each particle we find the distance to the 20th nearest 
particle and then calculate the overdensity inside a sphere of this radius. If the resulting fluctuation 
is above a certain critical value (in this case 0.5) we add the particle to the list of galaxies. This 
has been done for the F = 0.5 CDM simulation, for outputs with linear theory mass erg = 0.5 and 
as = 0.67. As the prescriptions are based on the overdensity of a given mass of particles, we label 
the resulting results “mbias 1” and “mbias 2” respectively. For the “mbias 2” prescription, we also 
add a random 30% of “field” particles to the list of galaxies, to represent peaks which could have 
formed in lower density regions but did not due to the finite resolution of the simulation. The 
second prescription involves finding the overdensity in a fixed volume about each particle (in this 
case a sphere of radius 2 /i“^Mpc). We add the particle to the list of galaxies if the overdensity 
is above 1.0 for one prescription or 3.0 for another. We label these sets “vbias 1” and “vbias 2” 
respectively. The F = 0.5 CDM simulation with as = 0.5 was used. 





Fig. 4.— UpkTr for the 4 different galaxy biasing prescriptions (see Section 2.3 for details) applied 
to simulations of the Standard (F = 0.5) CDM scenario. 

In all these cases, galaxy particles are preferentially selected in high density regions and 
are more clustered than the underlying mass. In Figure (^) we show how the bias factor 

galaxies /(y{r) mass Varies as a function of scale. We can see that all the models have some weak 
scale dependence on small scales, with the mbias 1 galaxies having the largest bias factor on all 
scales. 

The effect of these biasing prescriptions on the peak space density can be seen in Figure 4, 
where we also plot the linear theory prediction for the F = 0.5 CDM model. For filter scales 
above ~ 6 /i“^Mpc, the results are in good agreement with linear theory, as was the case for the 
unbiased simulations. Below this, we see a dehcit in the number of peaks, which is greatest for 
the mbias 1 and mbias 2 prescriptions, in which galaxies are selected according to the overdensity 
of a given number of surrounding particles and not governed by a physical length scale. For all 
biasing prescriptions, the deficit of peaks on small scales is something we should expect, as when 
we choose galaxies, we are effectively picking out peaks above a certain height in the evolved mass 
distribution. The peaks that fall below the threshold will disappear. The effect does not seem to 
be too drastic though, and on all but the smallest scales is smaller than the differences between 
cosmological models (see Figure 2). 
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3. Reconstruction of the power spectrum slope 

All the models we have tested in Section 2 started from initial conditions with Gaussian 
random phases. In these cases, all statistical information about the initial density is contained 
in the power spectrum, P{k). From Equations 1 and 2, we can see that the shape (but not 
the amplitude) of the linear power spectrum is responsible for the space density of peaks in the 
initial density field. In this Section we will present a simple and stable method for inverting these 
relations and recovering the shape of P{k) from the space density of peaks. We will then apply 
this method to the space density of peaks measured from simulations. Since Section 2 has shown 
us that on quasilinear scales this is the same as Upk in the initial conditions, we will therefore be 
reconstructing the initial power spectrum shape. 

Given the number density of peaks for a Gaussian field in equation (1), what can we say about 
its power spectrum? In other words, what information about P{h) can be obtained from the ratio 
of its moments, i?*, in equations 2-3? It is clear from these equations that i?* is independent of 
the amplitude of P{k), but it depends on its shape. If we wanted to compare with a known family 
of power spectrum shapes (e.g. CDM), we could just find the best parameter in the family (e.g. 
r) to fit the values of i?* as a function of rj. To hnd the power spectrum shape more generally, 
one could try using a numerical inversion technique (eg. a generalization of the method of Lucy 
1974) to extract the shape of P{h) from the measured values of i?*. 

A simpler approach, which we adopt here, is to approximate P{k) locally with a power law 
P{k) = k^: 


_ dlogP{k) 

If the slope n{k) is constant with scale, then we substitute P{k) = /c” exp [—{kRf)'^] into 
Equations 2 and 3 and find (Bardeen et ah 1986): 

7 6 \i/2 

This is only exactly true if n is a constant. However, if we assume that P{k) has a slowly 
varying shape, then it will remain a good approximation when applied locally, n is then given 
by the local value at k ~ l/r/, (where the Gaussian filter gives its maximum contribution). In 
Appendix A we show how to find the best effective relation between k and Vf. We can then use 
Equation 4 to translate the estimated values of R* at a given r/ to mean values of n at a given 
k. In Appendix A, we also test the method with various different analytic power spectra, going 
from the power spectrum to the peak density and then back to the reconstructed power spectrum 
slope. We hnd that the method is stable and works well. 
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3.1. Reconstruction from the peak density in simulations 

When a density field undergoes non-linear evolution under gravity, mode coupling causes a 
relative transfer of power between large and small scales. In the cases considered here, this has 
the effect of making P{k) less steep on small scales (see, eg. Baugh &: Efstathiou 1994). In Figure 
9 we show n{k), the logarithmic slope of P{k) measured from the evolved mass distribution of our 
three different cosmological models. This was estimated by first calculating P{k) using an FFT 
and then applying a two-point finite difference operator to recover n at the points shown in the 
figure. We also plot the linear theory predictions for n, which enables us to see that non-linear 
evolution of n is visible in the simulations on scales k > 0.07 /i“^Mpc. From Equation Al, using 
r] ~ 0.4, this scale is equivalent to r/ < 20 /i“^Mpc. 

We now apply our reconstruction method described in to the peak density measured from the 
simulations. This results in an estimate of n as a function of scale which we plot in Figure 10. 
The reconstructed n values are reasonably close to the linear theory values, as we would expect. 
Again, the T = 0.5 CDM simulation suffers the most from non-linear effects. It is interesting that 
the non-linearities (ie. merging) which affect the peak density cause it to be lowered, which means 
that the recovered n is slightly steeper than linear theory. We can see this on the smallest scales 
in Figure 11, in which we plot the reconstructed n for different output times in the T = 0.5 CDM 
simulation. This is opposite to the effect of non-linear evolution on the directly measured value of 
n (see Figure 9). 

In Figure 12 we show that the reconstruction also works well on the peak density measured 
from the biased simulations, and again works best on larger scales. 


4. Discussion and conclnsions 

Although a local transformation of a density field can change the heights of peaks, it should 
not affect the peak number density. Thus, as far as the effects of both gravity and biasing are 
approximately local one should expect the galaxy peak density to be a preserved quantity. Gravity 
is local on linear scales and there has not been enough cosmic time for gravity to affect (or become 
non-linear on) scales much larger than ~ 8/i“^Mpc. For similar reasons it also seems difficult for 
biasing to be significantly non-local (see also Gaztahaga & Frieman 1994 and references therein). 
Of course, it may also be possible to have non-local transformations which preserve the peak 
number density. 

Even without attempting to reconstruct the power spectrum shape from it, the space density 
of peaks appears to be an interesting statistic. At the present time, the two cosmological models 
which are the most favoured alternatives to T = 0.5 CDM are Low density CDM and MDM, which 
we have seen have very differently shaped linear power spectra and different peak space densities. 
The similarity of their shapes on quasilinear scales, and the uncertainties from biasing on smaller 
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Fig. 5.— The measured non-linear power spectrum logarithmic slope, n for different cosmological 
models. The points are the mean of results for 5 simulations, with the error on the mean calculated 
from their standard deviation. The curves are the linear theory values for n. 
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Fig. 6.— Power spectrum logarithmic slope , n reconstructed from the peak space density measured 
from simulations of three different cosmological models. 
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Fig. 7.— Reconstructed power spectrum logarithmic slope, n for three output times of the F = 0.5 
CDM simulation. 



k (h Mpc^^) 

Fig. 8.— Reconstructed power spectrum logarithmic slope, n for the different galaxy biasing 
prescriptions descibed in Section 2.3 and applied to the F = 0.5 CDM simulations. 
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scales mean that galaxy clustering observations generally favour neither one or the other. We 
propose that observational measurements of the peak space density would constitute a simple 
means of telling them apart. 

Attempts have already been made to constrain models based on linear power spectra 
reconstructed from the non-linear observations using mapping formulae calibrated against Ai-body 
simulations (as mentioned in Section 1). Unfortunately, two different analyses have come up with 
different conclusions, each favouring a different model. Low density CDM for Peacock & Dodds 
(1994) and MDM for Baugh & Gaztanaga (1996). This may be because of the sensitivity of the 
non-linear mapping formulae to cosmological parameters such as D, or the fact that they work 
less well with steep power spectra. The reconstruction from the peak space density should not be 
sensitive to the values of D and A, for example, and should be more insensitive to variations in 
galaxy biasing. It should be borne in mind, though, that the work in this paper is aimed at the 
quasilinear regime and that the non-linear mapping formulae can, in principle, be used as a tool 
to study much smaller scales (although the uncertainties due to biasing could be large). 

If one assumes that the initial conditions were Gaussian, we have shown that it is possible to 
reconstruct the shape of P(k) on interesting scales from the peak density. The method appears to 
work best for the models with the steeper slopes on small scales, MDM and P = 0.2 CDM. This 
is just as well, as these models are more favoured by galaxy clustering data. The effects of sparse 
sampling, boundary conditions and redshift distortions must be studied in conjunction with an 
application to real data. As with the genus statistic, a large contiguous volume is necessary in 
order to estimate Upk reliably. With the next generation of large redshift surveys (Colless 1997, 
Gunn & Weinberg 1995) it should be possible to measure Upk with comparable accuracy to the 
simulations we have presented here, and probably with even better accuracy on large scales. 

In conclusion, we have shown that the space density of peaks in a density field smoothed 
with a Gaussian filter is independent of time, following the linear theory prediction, even on scales 
where the clustering is mildly non-linear. We have also shown that the peak space density in 
itself can be used to distinguish between cosmological models, including the currently popular 
MDM and Low density GDM scenarios. The peak density could also be used to compare with 
non-Gaussian models, if predictions become available. We have developed a simple, stable method 
of recovering the power spectrum from the space density of peaks and demonstrated that the 
method works using analytic power spectra and the simulation results. An application of the 
method to a contiguous, densely sampled galaxy survey should yield a useful quantity, the shape 
of the linear power spectrum on scales k ^ 0.4 /i“^Mpc“^. 
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A. Prom the peak density to the power spectrum slope 

We will try a simple local reconstruction method which assumes that the shape of P{k) changes 
slowly with scale. In this case, we can approximate P{k) locally with a power law P{k) = k'^. For 
the case of interest, i.e. a density field smoothed with a Gaussian filter the relation between i?* 
and rf for a power law P{k) is given by equation 4. We have seen in Section 3 that this is only true 
for all rj if n is a constant. If the shape is slowly varying, a good estimate of n will be the mean 
local value taken assuming a one to one mapping between values of k and rj. We need to find 
a good effective relation between k and rj in order to use equation 4 to translate the estimated 
values of i?* at a given rj to mean values of n at a given k. The form we use is as follows, 

k = r]—, (Al) 

where in general tj may have some dependence on the scale k which will be unimportant if the 
shape of P{k) varies slowly. For a given family of P{k) the value of r/ = 'r]{k) can be found by 
equating the value of the local slope n at k with the mean value n obtained from equation 4, after 
finding the exact relation for = R^,{rf). In the case of a simple two power-law model with a 
break: 


P{k) = Ak^ [1 + {kc/kr] , (A2) 

we find that rj is indeed constant as a function of k: 

= (A3) 

and only depends weakly in the shape (but not the position) of the break (at kc). 


In our reconstructions we take r/ to be a constant. We find the best value of rj for each given P{k) 
by finding the value which gives the smallest difference between the measured value of Upk and the 
value predicted from the reconstructed P{k). We test this for a few arbitrary power spectra which 
we plot in Figure The F CDM models are taken from Efstathiou, Bond & White (1992). 


In Figure 0 we show, as a function of rj, the rms percentage error per logarithmic bin in the k 
range corresponding to Rf = 2 — 80/i“^Mpc for; a) the reconstructed value of (n -|- 5) (continuous 
line), compared to the input slope, and b) the value of the number density Upk calculated from the 
reconstructed P{k) (dashed-line), compared to the value of Upk calculated directly from P{k). 


As expected, the two power-law model (|A2D gives zero error at the rj predicted by (^). For the 
other models the error does not approach zero, indicating that the true t] is not exactly a constant. 
Nevertheless the minimum rms error value (which defines the best value of rj) is typically small. 
When applying the method to observations, as the true P{k) is not available, we can only estimate 
T] by comparing the errors between the observed and predicted Upk (dashed lines in Figure 0D, 
but as shown in the Figure this recovers well the best value of r/ at the minimum rms. 
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Fig. 9.— Different shapes for P{k): CDM like models with F = 0.05 (continuous left line) and 
F = 5.0 (continuous right line), 2 power-law model (|A2|) with n = —3 and m = —4. 



Fig. 10.— For each of the models in Figure 7 we show the percentage rms error per logarithmic 
bin in: a) {n + 5), where n is the logarithmic slope of P{k) (continuous line) b) the number density 
of peaks (dashed line). 
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In Figure we show a comparison of the reconstructed logarithmic slope i) for the values of rj 
that minimise the differences between the input and reconstructed Upk- The k range corresponds 
to Rf = 2 — 80/i“^Mpc, which is the input range in Upk- The agreement in all cases is quite good 
(less than 0.5%rms error per bin). 

In Figure we show two ad-hoc example power spectra where the method does not work quite 
as well as the changes in slope with scale are more rapid. The models shown correspond to power 
spectra with two sharp breaks or bumps. The rms error in the reconstruction is larger (about 
1.5%rms per bin), but note that the main features, such as the positions and relative amplitudes 
of the breaks, are still recovered. In these extreme cases a more careful reconstruction method 
should be used. 


Once we have the shape n = n{k), we can of course reconstruct P{k) itself up to an arbitrary 
constant by carrying out a simple numerical integration. This can be seen in Figure 13 where 
we show a comparison between the reconstructed and the input P{k) for the values of rj that 
minimized the differences between the input and reconstructed Upk- To reconstruct P{k) the 
amplitude has been matched so that cr^Rf = 8) = 1 with Gaussian smoothing (we have also 
rescaled some of the lines in the figure for clarity). As can be seen in this figure, the agreement in 
all these cases is quite good so that P{k) can be easily recovered, up to a constant. 
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Fig. 11.— The reconstructed (dashed-line) and the original (continuous line) slope of the power 
spectrum for different CDM like models with F = 5.0 (top), F = 0.5 and F = 0.05. The 2 power-law 
model ([A2|) with n = —3 and m = —4 is also shown. 



Fig. 12.— The reconstructed (dashed-line) and the original (continuous line) slope of the power 
spectrum for two models with rapid slope variation. 
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Fig. 13.— For each of the models in Figure ^ we show the reconstructed (dashed-line) and the 
original (continuous line) power spectrum. 




